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Superconducting gap structure of the 115's revisited 
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Density functional theory calculations of the electronic structure of Ce- and Pu-based heavy 
fermion superconductors in the so-called 115 family are performed. The gap equation is used to 
consider which superconducting order parameters are most favorable assuming a pairing interac- 
tion that is peaked at (7r,7r,5z) - the wavevector for the antiferromagnetic ordering found in close 
proximity. In addition to the commonly accepted d^2_yi order parameter, there is evidence that an 
extended s-wave order parameter with nodes is also plausible. We discuss whether these results are 
consistent with current observations and possible measurements that could help distinguish between 
these scenarios. 

PACS numbers: 74.20.Mn, 71.27.-(-a, 74.70. Tx 



^^ CeCoIns is a heavy fermion superconductor [TI , which has been shown to lie in close proximity to an antiferromagnetic 

OQ quantum critical point PHI]. Its structure consists of layers of square planar Ce atoms, similar in that respect to 

^_^ the cuprate superconductors. Soon after the discovery of superconductivity in CeCoIus many measurements were 

(—| performed that were consistent with lines of nodes in the superconducting gap, including specific heat, thermal 

conductivity [5], spin-lattice relaxation rate [6^, and penetration depth [7 . Ideally, ARPES measurements could 
^ reveal the anisotropic gap structure in fc-space, but current equipment has neither the energy resolution nor the 
5-H temperature range to perform such a study. In its absence probes which are directionally weighted averages over 
^-1 the Fermi surface have been used to identify the location of the nodes, including field-angle dependent specific heat 
^ and thermal conductivity [HHII], ffc2(^)[Il] point contact Andreev reflection measurementsflB] . and vortex lattice 

1 t structure [Ml [15] all of which now argue that nodes lie along the F — >■ (7r,7r,gz) direction in the tetragonal Brillouin 

2 zone. In addition, a neutron scattering resonance has been observed at (7r,7r,7r) which can be interpreted as consistent 
^ with d^2_y2 symmetry [IS]. Taken together the experimental evidence appears fairly compelling that CeCoIus is a 

'■ superconductor with a d^^-yi gap structure. Measurements on other members in this family are relatively scarce 

^-H due to the fact that either pressure is required for superconductivity to occur (as in CeRhlns, Celna, Ce2RhIn8, and 

Q CePt2ln7[T7H2n]) or there are radioactive considerations (for PuCoGas, PuRhGas, and PuCoIus [5TH1S] ) . However, 

O with the exception of a few measurements '24! those measurements which have been made on these compounds are 

identical to the results found for CeCoIns 25-30 . 

,-_l On the theoretical side, the BCS theory of superconductivity allows one to calculate the superconducting gap 

J> structure at T = K by solving the gap equation: 

^ A(k).-VF(k,k')^^E^, 



^1^ where A(k) is the superconducting gap function, ek is the electronic dispersion, and F(k, k') is the pairing interaction. 

Z^ To solve this equation the low energy electronic structure (ie. the Fermi surface, and effective masses) and the 

,__! pairing interaction F(k, k') must be known. This would appear to be a formidable if not impossible task for such 

|L* strongly correlated materials such as the cuprates or the heavy fermion materials which we wish to studv[5Tj. It is 

. ^H well known that for strongly correlated materials there are no well controlled approximations for calculating their 

^^ properties. However, it does appear from several cases that simple density functional theory (DFT) calculations 

J-j provide a reasonable starting point for computing Fermi surfaces in good agreement with experiment, even if the 
mass renormalization can be off by orders of magnitude |32H34j . The computation of the pairing interaction is more 



problematic. In weak coupling spin fluctuation theory, however, one may approximate r(k, k') by a renormalized 
spin susceptibility C/x5(q, a;)C/ where U is the screened Coulomb potential [55H55] . Whether or not this is rigorously 
true in the strongly correlated materials of interest, such as the cuprates, pnictides, and heavy fermions, remains 
a subject of debate, but at least in the case of the cuprates and the iron based superconductors, the observed gap 
phenomenology can be reasonably described utilizing this assumption [40l [41] . 

More specifically, from either a strong coupling, weak coupling, or experimental perspective the spin susceptibility 
of the cuprates is peaked at {tt,tt). Given the LDA determined Fermi surface with the largest density of states near 
(±7r,0) and (0,±7r) the order parameter which will naturally maximize the superconducting gap is a d^2_y2 order 
parameter. In the Fe-based superconductors such as doped LaFeAsO and BaFe2As2 the spin susceptibility is peaked 



at (7r,0). In this case, the combination of a hole pocket centered at the T point, and an electron pocket at the X 
point (of the unfolded Brillouin zone) leads to the suggestion of the so-called s± gap structure [42l - [45] . Furthermore, 
weak coupling approaches can calculate the variations in the spin susceptibility for different dopings and crystal 
structures. Inputting the renormalized spin susceptibility and electronic structure leads to strong gap variations and 
often accidental nodes which are consistent with many of the experimental observations |41j 

Since most of the experimental evidence for d-wave symmetry in the 115 materials is indirect, the story of the 
Fe-based superconductors compels us to reconsider the possible superconducting order parameters in the 115 family 
of superconductors, which are also compensated multiband metals like the Fe-based materials. The goal of this 
work is not to conclusively show which order parameter is correct in the various 115 materials. Rather, we assume 
a BCS-like formalism for the superconductivity with a repulsive interaction that is peaked in the vicinity of the 
nearby long range magnetic order and explore whether the d-wave order parameters are the only likely candidates for 
superconductivity. We find that alternative scenarios are at least plausible particularly for the Pu-based materials. 
For completeness it must be stated that it is also not obvious whether a BCS picture such as this is even the correct 
starting point to understand superconductivity in these materials, or whether other q vectors should be considered 
for the pairing interaction. However, in the absence of a microscopic theory the phenomenology at least is consistent 
with what is known about cuprate and Fe-based superconductors, and if it can be shown to be consistent in the heavy 
fermion superconductors as well, it provides a methodology by which one can tailor higher T^s within a family of 
superconductors without detailed knowledge of the pairing mechanism. 
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FIG. 1: (color online) (top) Three Fermi surface sheets of CeCoIns. (bottom right) The Fermi surfaces of PuCoIns. 



The WIEN2K code|46| is used to calculate the electronic structure using the generalized gradient approximation 
and adopting the Perdew-Burke-Ernzerhof exchange correlation potentials I47J. We included spin orbit interactions 
on the /-electrons with a second variational method. 

In the Ce-based 115s the evidence is fairly clear that the spin susceptibility is peaked at {Tr,Tr,qz). For Celns, 
CeRhlns, CeCoIns, Ce2RhIn8, and CePt2ln7, it has been shown that a dome of superconductivity emerges around 
a quantum critical point associated with the suppression of an antiferromagnetically ordered phase with Qaf = 
f7r.7r.J) [48U54j . We expect a similar situation to exist for Ce2PdIn8 and Ce9 Coins [551 [56] . Furthermore, the neutron 
resonance observed in the superconducting state of CeCoIns [IB], is either a result of the gapping of the spin exci- 
tations present at Q — (7r,7r,7r)[S7] or a consequence of the pairing interaction being peaked at this wavevector[55]. 
Consequently, assuming a uniform U the pairing interaction is indeed peaked at Q — (77,71,(7^) for all of the Ce-based 
superconductors listed above. 

In Fig. [IJ we show the Fermi surface for CeCoIns. It is a very complicated 3D Fermi surface. Tight binding models 
often approximate the Fermi surface in the 115s by multiple cylindrical electron Fermi surfaces centered at {'K^ir^kz) 
and a hole like Fermi surface centered at r(e.g. |59)). For the purpose of discussion a cartoon of the resulting Fermi 
surface is shown in Fig. 2 for kz—Q. 




(71,0) 



FIG. 2: (color online) Cartoon Fermi surface of the 115's in the {k^, ky) plane. The wavevector Q = (7r,7r) represents where 
the pairing fluctuations are likely maximal assuming a constant U and orbital matrix elements. Depending on which states 
they dominantly pair (ie. Qi or Q2) will dictate which order parameter symmetry is most favorable. 



We now ask the question of what gap symmetry is favored by such a Fermi surface with a pairing potential that is 
peaked at {n^n^qz). From Eq. (1) the order parameter will be chosen so that the maximal number of states within the 
gap energy of the Fermi energy which are separated by the wave vector Q = {n,TT,qz) can have an order parameter 
with opposite sign (ie. A(k) = - A(k + Q)). As in the cuprates, this results in a dx2_y2 symmetry when a large DOS 
exists at (7r,0,fcz) and equivalent points. This certainly is a possibility for the 115s as emphasized in Fig. 2 by the 
vector Qi. However, a second possibility is that the fluctuations strongly pair the states on different Fermi surface 
sheets (e.g. those sheets labeled a and e in Fig. 2 with vector Q2). This situation is analogous to the Fe-based 
superconductors [12H1S]- Consequently, an s± type pairing symmetry, belonging to the Aig representation where the 
order parameter changes sign between different Fermi surface sheets, would be favored. Given the complexity of 
the actual Fermi surface, and the need to overcome the strong onsite Coulomb repulsion U [501 ISI] it is likely that 
accidental nodes exist if the s± order parameter was the dominant order parameter. If the coupling between Fermi 
surface sheets is sufficiently weak then it is also possible to have different symmetry representations on different Fermi 
surface sheets (ie. s± on a and 7 and dx2-y2 on /3) [62ll63j. 

As an initial investigation to determine whether an s± or a d^2_y2 order parameter is more likely to occur in any 
particular 115 compound we examine the band structure more closely. One should consider the states which once 
renormalized to account for the enhanced Sommerfeld coefficient of the specific heat are within the spin fluctuation 
energy of Ep. dHvA measurements show that the experimentally measured Fermi surfaces are in good agreement with 
those computed using non-interacting DFT theory [32) . Thus, even though the mass renormalization is not correct, 
the DFT computed electronic structure should give a reasonable first guess as to the leading pairing instability. We 
suggest that the dx2_y2 order parameter will be dominant when the majority of states are found near (7r,0,/cz) and 
equivalent points, while the s± order parameter will be favored when the majority of states are found near the (0,0, /cz) 
and (7r,7r,fcz) points. For completeness all representations consistent with the tetragonal symmetry of the lattice should 
be considered |64j, but given the crudeness of this analysis for simplicity we only consider the d^^-y^ and s± order 
parameters as the most likely scenarios. 

Fig. 3 shows the band structure for several Ce-based superconducting materials, which are consistent with other 
band structure calculations where available [3H [SS]. For the 115s there are states near F and a large cylindrical 
Fermi surface centered at (7r,7r,fcz) which could favor an s± order parameter, but especially when considering the 
entire kz dispersion, there are significantly more states available for superconductivity near (7r,0,fcz) and equivalent 
points which would indeed favor dx2-y2 pairing. A non-uniform renormalization in fc-space or strong variations in 
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FIG. 3: (color online) (a) Band structure from top to bottom of CeCoIns 
points reflects the amount of /-character, (b) Zoom of panel (a). 



Ce2CoIn8, Ce2PdIn8 and CePt2ln7. The size of the 



orbital character could easily shift the preference to an alternative order parameter. dHvA results indeed indicate 
that the mass renormalization is nonuniform across momentum space. However, they show that the small pocket 
around the F point is actually the lightest band [32 which thus further supports the d^^-y^ gap structure over s±. 
While theoretical studies of the order parameter in the 115's generally favor (1^2 _y2 symmetry [SHIIMIj a study of the 
orbitally degenerate Hubbard model with tight binding parameters for the 115s also found alternative representations 
depending on the Coulomb and crystal field parameters WT . 

We also show the band structures for Ce2PdIn8, Ce2CoIn8, and CePt2ln7 in Fig. 3. The 218's can be visualized 
as a bilayer variant of the 115 structure. Indeed near 0.2 eV a doubling of the number of bands is clearly visible 
(especially clear near (TTjTTjfc^)). However, the number of Fermi surface sheets remains unchanged and the electronic 
structure near Ep is qualitatively similar between Ce2 Going and CeCoIns. The electronic structure of Ce2PdIn8 near 
Ep is similar to Ce2 Going with a small shift of the chemical potential reflecting the extra electron of Pd relative 
to Co, although it is clearly not a rigid shift. Finally, for CePt2ln7, the additional Ptln2 bands relative to the 115 
structure does result in additional Fermi surface sheets, although quite similar in structure to the other large Fermi 
surface sheets [BB]. While the details vary, the main features observed in the band structure of CeCoIn5 persist in 
these three other compounds as well. Namely, the dominant weight is found close to the {-Kfi^kz) points suggesting 
dominantly dx2_y2 gap symmetry. The magnetism and superconductivity present in the 115's appears robust to these 
minor changes in the electronic structure. 

In Fig. 4 we show the analogous band structures for the Pu-based superconductors, again in agreement with 
previous publications |55H7^ . As seen in Fig. 1, the Fermi surface of PuCoIns is qualitatively similar to CeCoIns, 
with respect to cylindrical- like Fermi surface sheets centered at (7r,7r,(72), a small hole pocket near F and a large 
Fermi surface sheet. However, one can also tell that in PuCoIns, the hole pocket is larger, the cylindrical sheets are 
closer to {TT,T:,qz), and there is less weight near {T:,0,qz) and equivalent points. Note, that in the case of the Pu-based 
superconductors it is much less clear where to guess the pairing fluctuations might be peaked since the analogous long 
range order has yet to be found. However, for the purpose of this discussion we will again assume that the pairing 
fluctuations are peaked at {T:,TT,qz). This is indeed what one would assume on the basis of a uniform U spin fluctuation 
calculation [73]. In contrast to the case of the Ce-based superconductors the arguments for dj;2_y2 symmetry are much 
less conclusive. Relative to the Ce-based materials there are fewer states in the vicinity of {7T,Q,kz) and more states 
near the Fermi level at (7r,7r,fcz). This leads us to speculate that the potential for an s± order parameter is greater in 
the case of Pu-based 115s than in the case of Ce-based 115s. Note also that in contrast to the Ce-based compounds. 



the bare / level is situated directly at the Fermi level in these Pu compounds. As a consequence, the bilayer variants 
such as Pu2CoGa8 have dramatically different electronic structure relative to FuCoGas. This is a possible reason for 
why superconductivity is much less prevalent in the relatives of the Fu-based 115s [721 El]; while it is much more 
robust in the Ce-based 115's. 
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FIG. 4: (color online) (a) Band structure from top to bottom of FuCoGas, PuCoIus and Pu2CoGa8 (which is not supercon- 
ducting). The size of the points reflects the amount of /-character, (b) Zoom of panel (a) 



We now briefly reconsider the experimental evidence for the position of the nodes in these materials. Ideally, field 
angle specific heat and thermal conductivity measurements can give the position of the nodes. However, they rely on 
the ability to accurately model the true electronic structure, which as we have emphasized, is poorly understood in 
heavy fermion materials. Furthermore, a similar 4-fold modulation has also been reported in some of the Fe-based 
superconductors [75J [75] . Similarly, while the point contact Andreev spectroscopy data could differentiate in favor of 
a dx2-y2 order parameter relative to a dxy order parameter, more work is needed for a realistic electronic structure 



whether alternative gap structures including s± with accidental nodes are possible [131 El] ■ We also find it interesting 
that the universal limit of the residual term in the thermal conductivity |78j is not obeyed with La-doping in CeCoIus 
[75] . This likely refiects the multiband nature of superconductivity in these materials [50] . 

What experiments could help resolve the question of the order parameter in these systems? Based on similar work 
which is being proposed to differentiate various gap symmetries in the Fe-based superconductors it is easy to propose 
similar measurements for the Ce- and Fu-based 115 superconductors, even though we acknowledge that temperature, 
pressure, and radioactive constraints may make many of these measurements very difficult. Certainly, phase sensitive 
experiments as done on cuprates [811182) would clearly be the most definitive tests of the order parameter. Additional 
studies such as examining details of impurity states with scannning tunneling microscopy ^83, and the quasiparticle 
interference pattern[53|, the neutron resonance [5S], and Andreev bound states [55] are a few examples of probes 
that will further help distinguish between different possible order parameters. In the exceptional case of the 20 K 
FuCoGas superconductor, ARFES measurements are obviously also potentially enlightening. Note, that one should 
also consider the possibility that the order parameter may change (and could gain/lose nodes) as a function of tuning 
parameters such as pressure or doping [5T, as has been suggested in the Fe-based superconductors. 

We have examined the band structures for various Ce- and Fu-based superconductors within the 115 family. Based 
on the assumption of pairing fiuctuations which are peaked at q^ = Qy — tt we find that the d^^-y^ order parameter 
is the likely candidate for superconductivity in the Ce-based materials. However, we also find that the s± order 
parameter may be favored if an anisotropic mass renormalization favors the gapping between states in the vicinity 
of (0,0, q^) and (7r,7r,g^). Experimentally this does not appear to be the case in CeCoIns, but it is a much stronger 
possibility for the Fu-based superconductors. The purpose of this survey is to stimulate further discussion of the 
order parameter in the 115 based superconductors. Future theoretical and experimental work is needed to clarify the 
order parameter in these materials. 
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